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LONG-TERM  GOAL 

Presently,  all  US  ship  mine  vulnerability  studies  are  based  on  the  restrictive  assumption  that  the 
ship  is  traveling  through  a  calm  sea.  Our  goal  is  to  remove  this  restriction  and  thereby  allow  a  more 
accurate  assessment  of  ship  vulnerability  to  sea  mines. 

OBJECTIVES 

We  are  developing  a  model  to  predict  the  pressure  field  surrounding  a  ship  advancing  in  a 
seaway.  This  pressure  field  consists  of  a  steady  component  due  to  the  forward  motion  of  the  ship  and 
an  unsteady  component  due  to  the  oscillatory  motions  induced  by  the  incoming  waves.  We  are 
particularly  interested  in  predicting  the  pressure  signature  on  the  seafloor,  because  this  is  where 
pressure-sensing  mines  are  typically  located.  Since  the  ship  pressure  field  decays  with  depth,  high- 
frequency  components  of  that  field  are  filtered  to  a  greater  extent  as  depth  increases.  Additionally,  the 
mine  samples  the  ambient  pressure  field  at  varying  but  relatively  low  frequency.  Both  of  these 
observations  suggest  the  need  for  a  model  that  includes  a  seafloor  and  is  accurate  at  low  frequencies. 
Most  conventional  strip  theories  assume  the  fluid  to  be  infinitely  deep  and  are  valid  only  for  high 
frequencies.  They  are  not  appropriate  for  our  purpose,  which  is  why  we  have  chosen  a  three- 
dimensional  formulation. 

APPROACH 

Our  first  step  was  to  derive  the  finite-depth  3D  Green  function  corresponding  to  a  translating 
source  with  an  oscillating  strength.  We  confined  our  attention  to  simulations  in  the  frequency  domain. 
These  solutions  assume  a  sufficient  time  has  passed  since  the  ship  first  encountered  the  waves,  and  all 
motion  transients  have  dissipated.  The  frequency-domain  formulation  is  computationally  efficient  and 
therefore  more  suitable  for  vulnerability  assessment  models.  Next,  we  developed  an  algorithm  to 
numerically  evaluate  the  Green  function.  The  integrals  in  the  Green  function  were  evaluated  using  an 
adaptive  quadrature.  Before  we  could  implement  the  Green  function  into  a  panel  program,  we  needed 
to  prove  that  the  flow  field  can  be  represented  by  a  source  distribution  over  the  submerged  hull  surface 
and  along  a  contour  where  the  hull  intersects  the  free  surface.  This  representation  is  equivalent  to  that 
in  the  infinite-depth  case.  Once  the  flow  field  can  be  represented  by  a  source  distribution,  a  panel 
program  can  be  developed  to  compute  the  motions  of  a  ship  in  waves  and  the  resulting  pressure  field. 

WORK  COMPLETED 

Derived  the  finite-depth  3D  Green  function  for  a  translating  and  oscillating  source. 

Developed  a  numerical  algorithm  to  compute  the  Green  function  and  its  derivatives. 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

2.  REPORT  TYPE 

4.  TITLE  AND  SUBTITLE 

Prediction  of  the  Pressure  Signature  of  a  Ship  in  a  Seastate 


6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Surface  Warfare  Center, Coastal  Systems  Station, Code 
Rll, Panama  City, FL, 32407 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  also  ADM002252. 

14.  ABSTRACT 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

4 

3.  DATES  COVERED 

00-00-1998  to  00-00-1998 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Showed  that  the  flow  field  can  be  represented  by  a  source  distribution  over  the  submerged  hull 
surface  and  along  a  contour  where  the  hull  intersects  the  free  surface.  This  representation  is  the  same 
as  in  the  case  where  the  depth  is  infinite. 

Developed  a  panel  program  that  incorporates  the  Green  function  and  computes  the  motions  and 
pressure  field  of  a  ship  advancing  in  waves.  The  program  uses  the  above  representation  of  the  flow 
field. 


RESULTS 

The  finite-depth  3D  Green  function  at  the  field  point  (x,y,z)  corresponds  to  a  translating  source 
with  an  oscillating  strength  of  frequency  o  located  at  (^,r|,Q  and  is  expressed  as 
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where  pv  indicates  the  principal  value  of  the  integral  and 


h  =  fluid  depth 
R2  =  (x-q)2  +  (y-p)2 
r2  =  R2  +  (z-02 
r22  =  R2  +  (z+^+2h)2 

J0  =  Bessel  function  of  the  first  kind  of  order  zero 
U  =  speed  of  the  source 


The  wave  frequency  co  is  related  to  the  wave  number  k  through  the  dispersion  relation 

co 2  =  gk  tanh(k/c) 


where  g  is  the  gravitational  constant.  The  group  velocity  Cg(k)  is  also  a  function  of  k  and  is  defined  as 
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dk 


The  function  F(k)  is 


cocosh(k(z  +  /z))cosh(k(C,  +h))cos(k(y  -ri)sin0) 
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For  0  <  0  <  %/2,  the  poles  kx  and  k0  are  the  two  zeroes  of  the  denominator  of  the  fifth  term  of  G,  i.e.,  k, 
and  k2  are  the  roots  of  the  equation 


(kU  cos0  +o)2  -co2  =0 

Here,  we  assume  that  kl  <  k2.  Similarly,  k3  and  k4  are  the  two  zeroes  of  the  denominator  of  the  sixth 
term  of  G.  They  are  the  roots  of  the  equation 


(kU  cos0  -o)2  -co2  =0 


The  poles  kD,  n=l,2,3,4  are  functions  of  0,  and  k3  and  k4  exist  for  0  e  [0,7t/2].  Thus,  the  range  of 
integration  of  the  sixth  term  and  the  last  two  integrals  is  from  0  to  n/2.  The  poles  k4  and  k~,  exist  only 

for  0  e  [0  ,71/2],  where  0  is  defined  below.  The  fourth  and  fifth  terms  of  G  differ  only  in  the  range  of 
0  integration.  However,  the  fourth  term  contains  no  singularities  and  can  be  integrated  numerically 
without  any  difficulties  whereas  the  fifth  term  contains  two  singularities  k]  and  k2. 

The  definition  of  0  is  rather  complex.  Its  value  depends  on  U,  h  and  o.  When  h  and  o  are 
given,  we  can  define  two  values  k  and  U  where  k  is  the  solution  of  the  equation 

co  (k)-kCg(k)  =  o 


and  U  is  given  as 


U  =Cg(k) 

For  U  <U  ,  0  =0,  and  the  fourth  term  vanishes.  All  the  remaining  integrations  with  respect  to  0  are 
from  0  to  lit 2.  For  U  >U  ,  0  =cos r\U  tU  ). 

We  validated  the  Green  function  and  the  numerical  algorithm  used  in  its  evaluation  by 
computing  the  steady  flow  around  a  fully  submerged  Rankine  body.  The  Rankine  body  was  chosen  for 
its  simplicity.  It  can  be  represented  by  a  single  source-sink  pair  and  approximates  an  underwater 
vehicle  fairly  well.  The  Green  function  for  the  steady  flow  problem  is  a  special  case  of  our  more 
general  translating  and  oscillating  Green  function.  It  is  obtained  by  setting  the  frequency  in  our  Green 
function  to  zero.  The  figure  compares  our  results  with  those  of  Sahin  et  al.  (1994).  The  Rankine  body 
in  this  example  has  a  length  to  diameter  ratio  of  7  and  a  length  to  depth  ratio  of  2.25.  The  hull 
centerline  is  located  at  mid-depth  (a=  z/d=0.5),  and  the  Froude  number  is  0.7.  The  figure  shows  the 
axial  distribution  of  bottom  pressures  directly  below  the  body  (y=0).  Excellent  agreement  is  observed. 


x/L 


IMPACT/APPLICATIONS 

The  primary  purpose  of  this  research  was  to  improve  the  Navy’s  ability  to  assess  ship 
vulnerability  to  sea  mines.  However,  our  model  can  also  be  used  to  design  improved  pressure-sensing 
mines. 

RELATED  PROJECTS 

Under  a  Coastal  System  Station  Internal  Research  task.  Dr.  Thai  Nguyen  is  investigating  the 
effects  of  a  compliant,  muddy  seafloor  on  the  motion  of  ships  in  waves.  A  different  Green  function 
was  derived  and  computed  in  this  case  where  there  exists  a  thin,  denser  fluid  layer,  representing  the 
mud  layer,  below  a  lighter  upper  fluid  layer.  However,  the  same  panel  program  is  used  in  both  projects 
with  the  exception  of  the  subroutine  computing  the  Green  function. 
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